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Abstract: We study cascaded quadratic soliton compressors and address 
the physical mechanisms that limit the compression. A nonlocal model is 
derived, and the nonlocal response is shown to have an additional oscillatory 
component in the nonstationary regime when the group-velocity mismatch 
(GVM) is strong. This inhibits efficient compression. Raman-like pertur- 
bations from the cascaded nonlinearity, competing cubic nonlinearities, 
higher-order dispersion, and soliton energy may also limit compression, 
and through reaUstic numerical simulations we point out when each factor 
becomes important. We find that it is theoretically possible to reach the 
single-cycle regime by compressing high-energy fs pulses for wavelengths 
A = 1 .0 — 1 .3 /im in a /3-barium-borate crystal, and it requires that the sys- 
tem is in the stationary regime, where the phase mismatch is large enough 
to overcome the detrimental GVM effects. However, the simulations show 
that reaching single-cycle duration is ultimately inhibited by competing 
cubic nonlinearities as well as dispersive waves, that only show up when 
taking higher-order dispersion into account. 
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1. Introduction 

Compression of optical pulses can be achieved by first inducing a nonlinear phase shift on the 
pulse by self-phase modulation (SPM) from the cubic nonhnear response of, e.g., an optical 
fiber [1]. The phase shift creates a chirp across the pulse, which means that compression can 
subsequently be achieved in a dispersive material (such as a grating pair). In cubic soliton 
compressors both the SPM induced chirp and the compression is achieved in the same material 
[2]. The self-focusing cubic nonlinearity requires anomalous dispersion to compress the pulse. 
Due to collapse problems this limits the energy of the compressed pulse, and the requirement 
of anomalous dispersion restricts the accessible wavelength regime for soliton compression. 
Recent progress has shown that cascaded quadratic soliton compressors (CQSCs) may ef- 
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Fig. 1. Numerical simulations of soliton compression of a 200 fs FWHM pulse in a BBO 
crystal pumped at X\ = 1064 nm with a soliton number of A'gff = 8. (a) and (b) Temporal 
and spectral components of the FW | t/i | ^ = /i //in in the stationary regime {Ak = 50 mm^ ' , 
/in = 59 GW/cm^). The pulse is compressed to 6 fs at the optimal compression point [z = 
Zopti dashed line, cuts in (e) and (f)]. (c) and (d) In the nonstationary regime (Ak = 30 mm^' , 
/in = 29 GW/cm^), a 17 fs pulse compressed pulse with trailing oscillations is observed. 



fectively compress liigli-energy fs pulses down to ultra-short few-cycle pulses [3-11]. Here the 
nonlinear phase shift is induced due to phase-mismatched second-harmonic generation (SHG), 
which acts as a cascaded quadratic nonlinear process. The pump, or fundamental wave (FW), 
experiences an effective SPM from the cyclic energy transfer to the second harmonic (SH). The 
advantage is that the effective cubic nonlinearity induced by the cascading process can be made 
self-defocusing because it is controlled by the sign of the phase-mismatch parameter [12-15], 
and therefore normal dispersion can be used to compress the FW. 

As this compressor scheme exploits an effective self-defocusing cubic term from cascaded 
quadratic effects, the compressor will naturally be affected by the self-focusing cubic nonlin- 
earity inherent to any transparent material. This detrimental cubic nonlinearity must be counter- 
balanced and then exceeded to achieve compression [3,6, 1 1, 14]. Furthermore, the propagation 
in a bulk medium having self-focusing cubic nonlinear response can avoid collapse problems if 
the self-defocusing cascaded nonlinearity is strong enough [16]. Thus, this compressor works 
in a bulk configuration even with multi-mJ input pulse energies [14]. 

Experimentally, compression of mJ pulses was observed at near-infrared wavelengths of 800 
nm from 120 fs down to 30 fs when the compression was done externally (using either a prism 
pair or a near-lossless block of bulk calcite) [14] and when working as a soliton compressor 
from 120 fs to 45 fs [3] and from 35 fs to 20 fs [5]. In an important recent advance, spa- 
tially uniform compression using super-Gaussian flat-top beams was demonstrated also at 800 



nm [7]. At 1260 nm compression to 12 fs (3 optical cycles) was achieved [6], while at telecom 
wavelengths compression down to 35 fs was observed [4, 8]. Numerical simulations have pre- 
dicted sub-2-cycle pulses, but it was clear already in the first experiments that group-velocity 
mismatch (GVM) was a limiting factor for compression [6, 14]. It was observed that in the 
so-called stationary regime clean compression is possible, while in the so-called nonstationary 
regime GVM distorts the compressed pulse too much to be of any practical use, and severe 
reductions in compression capabilities is observed. As an example of this, the numerical sim- 
ulations in Figs.[na)-(b) and Figs.[Tlc)-(d) compare the pulse compression performance under 
equal conditions in the two regimes. In the stationary regime a 6 fs compressed pulse is ob- 
served while the nonstationary regime the GVM effects are much stronger, resulting in a 17 fs 
compressed pulse with trailing oscillations. 

Significant progress in understanding this was recently made by using nonlocal theory (see 
Ref. [17] for a review on optically nonlocal media). The GVM-induced Raman-like term found 
previously [6, 18] was shown to originate from a temporally nonlocal response function [10]. 
The nonlocal behaviour appears when approximating the phase-mismatched dispersive SHG 
process in the cascading limit as a nonlinear convolution between the FW and a nonlocal re- 
sponse function [19]. An accurate prediction of when the system is stationary or nonstationary 
was presented, and the nonlocal theory predicted that an oscillatory chirp was created on the 
FW for short enough pulses in the nonstationary regime, which effectively limits the amount 
of compression achievable and qualitatively explained the trailing oscillations observed. It was 
also argued that the temporal time scales of the nonlocal response function had an influence on 
the final compressed pulse duration, but a systematic investigation was not made. 

On the other hand, a recent study showed that the performance of the CQSC can conve- 
niently be described by scaling laws involving an effective soliton number A^^^f = A^Ihg "-^Kerr' 
appearing as the difference between the quadratic and the cubic soliton numbers [11]. Since A'eff 
depends only on input parameters (such as parameters describing the input pulse intensity and 
duration and the material), the compressed pulse properties can be predicted using these scal- 
ing laws. Appropriate input parameters can then be found giving compression to single-cycle 
duration. However, neither the experiments nor the numerical simulations have ever observed 
single-cycle compression, so obviously there are higher order effects that prevent compression 
to reach such levels. In particular, both experiments and simulations found optimal phase mis- 
match values, where the best compression is observed. It is the purpose of the present theoretical 
and numerical analysis to discover what defines these optimal operation points. 

2. Propagation equations 

The SHG propagation equations in the slowly-evolving wave-approximation (SEWA) are used 
to study pulses in a bulk quadratic nonlinear crystal with single-cycle temporal resolution. Ne- 
glecting diffraction as well as the non-instantaneous cubic Raman response the dimensionless 
equations for the FW («i) and SH {oh = 2(0i) fields Uio{B, , t) are [1 1, 20] 

{id^ +b\)Ui + \M'\'I^NmGS\UlU2e'^'^ +nLJ[Ui {\Ui\^ + Bn\U2\^) = 0, (la) 
{id^ - id[2dr + D',,,s)U2 + \Ak'\'/^NsHGS'2Ufe-''"''^ 

+ 2n^N^J'2U2{\U2\"+Bfi-'\Ui\^)=0. (lb) 

Higher order dispersion (HOD) is included through the operator D'j = T!m=2''"^j'"^ 

with the dimensionless dispersion coefficients Sj'"' = k^l"\T^^^^\kf''^\m\)^^ and ^y'"' = 
d"'kj/dco'"\(o=o)j- Since kj = iijCOj/c is known analytically through the Sellmeier equations 

of [21], the exact dispersion Dj ~ kj{03) — (o) — fi)y)A:|'' — kj{o3j) is used in the numerics [11], 



corresponding to a dispersion order = iij is the refractive index, and the phase mismatch 
of the SHG process is AA: = A:2 — 2ki . The Kerr cross-phase modulation (XPM) term B = 2 for 
type SHG while for type I SHG B = 2/3 [11]. The time coordinate moves with the FW group 
velocity Vg.i = 1 /k[^\ giving the GVM term di2 | ~ 2- The equations are reported in 
dimensionless form, T = f/Tin, where Tin is the FW input pulse duration, £, = z/L^.i, where 
Ld 1 = T-^/\kf^\ is the FW dispersion length, and finally Ui = Ei/S'i„ and U2 = E^j \fhS\^. 
Here is the amplitude of the peak electric input field, « = «i/«2, <^'i2 ~ di2Tm/\k\ |, and 
Ak' = AA:Ld,i- This scaling gives the quadratic (SHG) and cubic (Kerr) soliton numbers [6, 11] 

_ ^D.l'g'inft^?'^eff ..2 „ ^D.l»Ken-.l<ginfiJl 

c^nin2\Ak\ c 

where t/gff is the effective quadratic nonlinearity, and WKen-j ~ 3Re(x'-'-')/8ny is the cubic (Kerr) 
nonlinear refractive index. A^shg might seem poorly defined in Eqs. ([T]i because of the factor 
|AA:'| '/^ in front of it, but the choice will become clear later Self-steepening is included through 
the operators S'j = 1 + i{(OjTin)^^ The combination of self-steepening and GVM implies 

that in the SEWA framework the SH dispersion effectively is given by D2 eff = ^2 + '^2^ I f 



2 dT 

led symbols in our notation i 

dimensionless form of the corresponding unprimed symbol 



where V = cd^2/('>2'^2\kf'^\ [11^20]. We stress that all primed symbols in our notation are the 



3. Nonlocal model: reduced equation in the cascading limit 

The CQSC works in the cascading limit, i.e., \Ak'\ ^ 1, and we henceforth take A^' > as to 
have a self-defocusing cascaded nonlinearity. In order to support soUtons we then need nor- 

(2) 

mal group- velocity dispersion (GVD), i.e., kj > 0. We now seek to reduce the full SEWA 
model ([T]i, in order to get some physical insight into the compression process. In Ref. [10] it 
was shown that in the cascading limit Eqs. ^ can be reduced to a single equation for the FW 



Ui + Ni,„Ui \Ui\'~ NiiicU^ I dsR± {s)Ut {^,Z-s)=0. (3) 



In this derivation both self-steepening and HOD were neglected (Sj = 1, = 2), but as shown 
later this can straightforwardly be relaxed. Additionally the Kerr XPM terms were neglected. 
This dimensionless generalized nonlinear Schrodinger equation (NLSE) shows that the cas- 
caded quadratic nonlinearity imposes a temporal nonlocal response on the FW, governed by 
the nonlocal response functions R± [examples are shown in Fig. |2ta,c)]. This model quanti- 
fied the previous qualitative definitions [14, 18] of the stationary and nonstationary regimes. 

9 (2) 

When > 2Akk2 the system is in the nonstationary regime and the oscillatory response 
function R must be used. The criterion can be interpreted using characteristic length scales, 
giving Lqvm < i'cohi'D,2/27r; when GVM dominates its length scale Lgvm = 7in/|c/i2| becomes 
shorter than one controlled by the product of the coherence length Lcoh — 7t/\Ak\ and the SH 

GVD length scale Ld.2 = 7in/l^2^'' I- When d^2 < 2Akk^'^ the system is in the stationary regime 
and one has to use the localized response function R^. We will now derive this result in details. 
In the cascading limit A^' ^ 1 the nonlocal approach takes the ansatz 

C/2(^,T)=(/)2(T)exp(-/Afc'^) (4) 



This ansatz is assuming that all the dynamics in the propagation direction of the SH is dom- 
inated by the phase mismatch, and the condition for making this ansatz is that the coherence 



length Lcoh = tt/AA: is much shorter than any other characteristic length scale. This is true in 
the cascading limit except when the FW is extremely short, in which case the GVM length 
Lgvm = Ti„/\di2\ can become on the order of Lcoh- Assuming A^Kerr^2 ^ ^'^^^Nshg we may 
discard the Kerr terms in Eq. ( fTbl i. and get an ordinary differential equation (ODE) 

5f ' ^ + id'J-^ - Kk'<H = Al/'/'NsHGUf (5) 

where for simplicity we have only considered up to 2nd order dispersion and ne- 
glected self-steepening and the SEWA correction to the dispersion. We will come 
back to this point later. Fourier transforming Eq. ^ and using the Fourier trans- 
form pair 02(ii) ■^i'PiK^) = (2;r)-^/2X!ldTe'"^02(T) and 02(t) = ^-^[02](t) = 
(2a:)-'/2£ldi2e-'"^02(fi) the ODE ^ becomes simply {Q?5?^ - D.d\2 + ^')hi^) = 
—Ak'^^^NsHG-^iUf] (n), so we can solve it in the frequency domain as 



-1/2 



5PQ?-d[2^ + Ak' 



(6) 



We now use the convolution theorem, so that in the time domain Eq. (|6]l becomes 



02(T) = -^/ dsRis)Uti^,z-s) (7) 



/Ak' 

where the nonlocal response function in the time domain is R{t) = .^^^[R[£l)\. This result 
shows that in the cascading limit the SH is affected by the FW in a nonlocal manner 
The next step is to calculate R{t). To do this it is convenient to rewrite R in Eq. (|6]l as 



R{D.) = {2k)-" \ ^''+'l>'^b (8) 



where we have introduced the dimensionless frequencies and the sign parameters 

i2;,=«';2/25f , £i;=|Afc752'')-n;?|'/2 (9a) 

=sgn[n;,] , Sh =sgn[Afc75f ' - D!^] (9b) 

The advantage of this notation is that it is now easy to characterize the roots of the polynomial 
in the denominator, which is important when we later will generalize to HOD . 

When if, = +1 Eq. (O becomes a Lorentzian centered in and with the FWHM 20!^, see 
Fig. |2b). The roots in the denominator of Eq. dH) are complex i2 = ± /£2'^. The temporal 
response function, R^(x), can readily be calculated by taking the inverse Fourier transforrrQ 

T 2 + 17,2 

= ^^^^2:^^'<^P(~""V'^«)exp(-|T|/Tf,) (10) 

where we have introduced the dimensionless characteristic nonlocal time scales 

T„ = |n;,ri = \25'^^/d'nl T, = ii;-' = |Afc75f -ii:2|-V2 (H) 

The localized nature of Eq. (fTOl i is shown in Fig.|2a), and Zj, controls the width of \R+\ while 
Xa is the period of the phase oscillations. We note that Eq. (|6]l is defined so /"^dT7?+(T) = 1. 



'Note: on dimensional form .'^(t) = R^t/Ti^^) /7i„, which is independent on Tj,, since one must replace T^.b with the 
dimensional equivalent t„ = Ta.bTin in Eqs. jl0ll2t . Instead in the frequency domain both R and are dimensionless. 
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Fig. 2. The nonlocal response functions in the (a,b) stationary regime (5^ = 1) and (c,d) 
nonstationary regime (.s/, = — 1) as calculated for the simulation parameters in Fig.[T] The 
cusp in the origin of \M±{t)\ is typical for a Lorentzian response function. The spectral 
content of Uf = sech^(f/riji) having 100, 20, and 10 fs FWHM duration is shown in (b,d). 




When Sfe = — 1, R{Q.) has two simple poles at fi = H', ± HJ, = Q!^, making R{Q.) diverge as 
shown in Fig.|2d). Then R-{t) ~ exists as a Cauchy principal value giving 

= exp(-w„T/T„)sin(|T|/Tfo) (12) 

In contrast to 7?+, this response function is not localized but instead is oscillatory as a conse- 
quence of the two poles in R (ii) [see the example shown in Fig.|2tc)]. 

Now, using Eq. ^ with the ansatz ^ and plugging it into Eq. ( [Tal l we arrive at Eq. (O un- 
der the aforementioned approximations. This essentially shows the central result of Ref. [10], 
namely that due to the cascaded quadratic nonlinearity the FW can be described by a general- 
ized NLSE with a nonlocal temporal response in the Kerr-like SPM term. In Ref. [18] the non- 
stationary regime was defined as when GVM effects prevent the build-up of a nonlinear phase 
shift, which when applied to soliton compression consequently results in poor compression [6]. 
Based on the nature of the nonlocal response functions, we can now clarify that the boundary 
to the stationary regime is when si, changes sign, which in dimensional units Kk = Afcsr, where 

Thus, to be in the stationary regime the phase-mismatch must be significantly large, AA: > l^k^x- 
When GVM is weak compared to the phase mismatch then = -f 1, and the response func- 
tion (/?+) is monotonously decaying in magnitude: the convolution in the Kerr-like SPM term 
in Eq. (O provides a finite temporal response. Therefore this must correspond to the station- 
ary regime. Instead when GVM is strong compared to the phase mismatch then i/, = — 1, and 
the response function (/?-) is oscillating and non-decaying: the temporal response from the 
convolution is no longer finite. Thus, this must correspond to the nonstationary regime. 



Finally, to include the influence of self-steepening in the nonlocal theory, Eq. iT% would 
have the operator S\ acting on all nonlinear terms and Eq. ^ would have the operator ^2 acting 
on the RHS. In frequency domain we may thus use a steepening-corrected response function 
R^^{Q.) = ^(i2)[l + (a)27in)^'ii]. This does not change the transition value AA:sr. Lastly, the 
NLS-like nonlocal equation ( fTTI i will have the operator S\ acting on all nonlinear terms. It 
should also be stressed that self-steepening can affect the Raman-like term in Eq. (T% [20], but 
this effect does not appear in the nonlocal model used here because this would require taking 
into account higher-order perturbation terms [i.e., making a more elaborate ansatz than Eq. ^]. 

4. The weakly nonlocal limit 

The nonlocal response in Eq. (O can be better understood in the weakly nonlocal limit, 
where the width of the nonlocal response function is much narrower than the width of Uf. 
The resulting simplified equation gives a better physical insight [22], and is important be- 
cause it governs the initial dynamics (until pulse compression makes Ui so short that the 
nonlocal response is no longer weak). We evaluate the convolution in the frequency domain 
/r„d5/?(i);7f(^ ,T-s)= J';^^dQ.e-'^^R{Q.).^[Uf] (ii) for convenience. In the weakly nonlo- 
cal limit R{Q.) is approximated by a 1st order expansion around Q. — 0, where ■^[Uf]{Q.) is 
non- vanishing. This holds when R{Q.) varies slowly compared to ^[U^]{Q.). In this case 



R{Q.)^[uf]{D.) ^. 



dR 

R{n = 0)+£l—\n=o 



(14) 



which in the time domain equivalently is R{s)Uf{^,T — s) c^R{s) [uf{^,T) —sdUf{^,T)/dT] . 
However, in the nonstationary regime the integration is done over two simple poles located on 
the H-axis. Using the residue theorem the frequency integral J'^^dQ.e^'^'^R{Q.),^[Uf]{Q.) can 
be evaluated as a contour integral, giving a contribution from the Cauchy principal value of the 
integral, and a contribution from deforming the integration contour around the poles on the real 
f2-axis. The residual contributions from the poles to the frequency integral are [23] 



p(T,t/i) ==/sgn(T) 



2 2T2Tfo 



(15) 



which consist of an oscillatory component in form of complex exponentials with frequencies 
each weighted by the spectral strength of Uf at that frequency. Thus, the influence of 
this contribution becomes important when the FW is short enough for its spectrum to cover 
the range, where Q.'_j_ are located, see Fig. |3d). Using Eq. ([Bll and (27r)'/2^(0) = 1 and 
(27r)'/2d^/dn|n=o = 2Q.'J{Q.'^ + Sh^f), the nonlocal convolution is 



(16) 



Now introducing the effective soliton number A^^^ = -^shg ^ -^Ken' © becomes 
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(17) 



The first term on the RHS is a GVM-induced Raman-like perturbation caused by the cascaded 
SHG nonlinearity. It is Raman-like due to the asymmetry of R± [10], stemming from the phase 
term exp(— WaT/ To) in Eqs. dTol ) and (fT2] |. It has the characteristic dimensionless time Tr^shg = 
4|£2;,|/(n;2 + st^'^) = l\d\^lN<! [6, 10,11, IS], wWch on dimensional form reads 



Tk.SHG = TR.SHGTin = 2|(ii2|/Afc 



(18) 



The direct dependence on the GVM-parameter di2 impHes that the Raman-like perturbation 
disappears in absence of GVM. Equation (lT\ was also derived for the stationary regime si, = l 
in Ref. [lojE obviously without the contribution f/[*p(T,[/i). 

Eq. ([TtI i is a very strong result: It states that in the weakly nonlocal hmit the effective soliton 
number A^eff can be used in the scaling laws of [11] to predict, e.g., the optimal compression 
point. Previously these scaling laws were thought to hold only in the stationary regime [10, 11]. 
The result also tells us that in the weakly nonlocal limit, a central observation of Ref. [6] is 
confirmed: for a given, fixed value of A^, the Raman-like effect of the first term of the RHS 
of Eq. ([T7] | becomes increasingly significant with increasing A^shG' ^^^^ limiting the possible 
compression ratio. However, it is now clear that in the nonstationary regime, the Raman-like 
distortion is accompanied by an oscillatory perturbation term Ulp{t,Ui) which also increases 
with A^sHG- ^^'■^ '■^^ stationary and nonstationary regimes the Raman-like distortions place a 
limitation on the maximum soliton order, but in the nonstationary regime both terms on the RHS 
of Eq. dTTI i distort the compression, and the combined effect is more severe, a result that is also 
confirmed by numerical simulations below. On the other hand, it is also obvious that for low 
soliton numbers these detrimental effects are weak and it is actually possible to observe clean 
compressed pulses even in the nonstationary regime, as observed by both numerical simulations 
as well as experiments [3,5]|1. This is also apparent in Eq. (ITsT l where the oscillations can be 
neglected if ^[C/f ] {^±) is vanishing, i.e., if the compressed pulse duration is long enough. 

In the stationary regime the characteristic time scale that decides the strength of the nonlocal 
response is tj, = ThTin- The weakly nonlocal approximation applies when f/, <C Af, where Af is 
the FW pulse duration. But when does it apply in the nonstationary regime? We know that the 
width Ai2 of J^[f/f] is Ai2 oc Af Referring to Fig. |2td) we must require that the positions 
of the two poles i2± be sufficiently far away from the frequency range, where J^[f/f] is non- 
vanishing, i.e., 3> Ai2. In physical units this implies that the weakly nonlocal limit in the 
nonstationary regime can be expressed by the requirement Af ^ tatb/\ta — 

Let us evaluate this requirement. It is important to notice from Eq. ( ITTI ) that diverges 
at the transition hk^^, see also Fig. Ha). Thus, in the nonstationary regime the requirement 
At ^ tatbl\ta — tb\ holds even for quite short pulses as long as f„ and ff, are not too similar This 
is generally true close to the transition AA:sr, while away from the transition f^ ~ f^ because Ak 
gets small, see Eq. dTTl i and Fig.[3la). In this case we can no longer be sure to be in the weakly 
nonlocal limit. In the stationary regime the system will initially be in the weakly nonlocal limit 
Af 3> tfj except close to the transition AA-'sr, where f/, diverges. It is interesting to notice that the 
weakly nonlocal limit in the stationary regime is well away from the transition Afcsr, while in 
the nonstationary regime it is close to the transition A^si - We finally remark that ta also may 
diverge when GVM is negligible. This implies that the factor e^"^"'^!'^" ~ 1, so R{'z) becomes 
real and symmetric. Thus, the 1st order correction on the RHS of Eq. ( [TT] ) disappears because 
the Raman-like perturbation vanishes (Tfl shg = 0) and a 2nd order correction must be made. 

5. Numerical results and discussion 

We will now apply the results of the nonlocal theory to understand realistic numerical simula- 
tions of Eq. (dJ. It is important to stress that the theory neglects Kerr XPM effects and that the 
coherence length is the shortest length scale in the system. This latter requirement implies that 
the system can initially be well described by the nonlocal theory, but as the pulse is compressed 
the GVM (and other length scales) can become so short that this is no longer true. Therefore the 
accuracy of the nonlocal model will not always be enough to accurately predict the outcome of 
the numerical simulations and the experiments. However, since the nonlocal model often will 



The factor on the RHS ofEq. (Tt) was unfortunately lost during the proofs in Eq. (12) of Ref. [10]. 
'These experiments were actually done in the nonstationary regime according to the nonlocal theory. 



be an adequate approximation for a large part of the propagation through the nonUnear medium, 
we can still use it to understand the temporal dynamics until that potentially happens. 
We will now show that there are two main categories of compression limitations. 

1. Effects limiting phase-mismatch range where compression is possible and efficient 

(a) In the nonstationary regime M. < Msr the oscillatory nonlocal response function 
implies that compression is inefficient unless the soliton order is very low. 

(b) Competing cubic nonlinearities pose an upper limit Akc [10, 11] beyond which 
A^eff < 1 always. Close to this limit detrimental cubic XPM effects are observed. 

2. Effects Umiting the compression for a selected phase-mismatch value 

(a) The effective soliton order A'eff = (^shg ^^Kerr)'''' controls in the weakly nonlocal 
limit the compression factor fc = Tin/Afopt = 4.7(A^eff ^ 0.86) [11]. 

(b) Nonlocal effects. In the stationary regime Afopt is limited by the strength of the 
nonlocal response function f^. In the nonstationary regime Afopt is limited by the 
characteristic time 7r shg of the GVM-induced Raman-Uke perturbation. 

(c) Propagation effects pertaining solely to the FW, such as higher-order dispersion, the 
Raman effect (negligible in nonlinear crystals) and cubic self-steepening. 

(d) Competing cubic nonlinearities necessitate large quadratic soliton orders A^shg- 
which increases detrimental nonlocal effects such as the Raman-like perturbation. 

The numerical simulations of Eqs. ([TJ were performed using a j3-barium-borate crystal 
(BBO) as the quadratic nonlinear medium. The phase mismatch was changed through angle- 
tuning of the crystal in a type I SHG configuration (implying B = 2/3, see [11] for further 
details), and we are interested in M > 0, for which GVD is normal and 12 < (so Sa < 0). 

Figure [3] summarizes simulations of pulse compression of a 200 fs FWHM pulse centered 
at Ai = 1064 nm: the FW pulse duration at the point of optimal compression Afopt (dark cir- 
cles) is plotted as the phase mismatch Ak is sweeped. The strength of the cascaded quadratic 
nonlinearity N^^^q °^ Ak^ ' , while the Kerr nonlinearity remains unchanged. However, in the 
plot we keep Nes ~ 8 fixed by adjusting the input intensity = ^«ieoCffi^0This implies that 
the scaling law [11] predicts equal compression everywhere (dotted orange line). While we do 
observe such a compression in the regime around A^ = 50 mm^' [this example is shown in 
Fig-[na)-(b)> and also in Fig.lD, away from this point the compression becomes sub-optimal. 

Before explaining the results in detail, we must mention that the nonlocal time scales plotted 
in Fig. |3la) are not those of Eq. (fTTT i; these only take into account up to second-order dispersion 
(dispersion order — 2), while in the numerical simulations the dispersion is calculated ex- 
actly from the Sellmeier equations [11], and subsequently corrected in the SEWA framework up 
to 30th order [11]. This poses a correction on the nonlocal time scales as well as the transition 
to the nonstationary regime, which was calculated numerically: we replaced the polynomial in 
the denominator of R in Eq. ^ with D2.eff (evaluated in frequency domain). The transition to 
the nonstationary regime happens when a root-pair switches from being each others complex 
conjugate to being purely real and nondegenerate. Then D.a,b and t^j, can be extracted from 
these roots. The transition to the stationary regime ( fTlT l is now simply found when f/, diverges. 
Fig.[3]shows for comparison also f^ calculated with ~ 2, i.e., using Eq. (fTTT i. 

The degrading compressor performance observed for large Ak in Fig. [3ja) is caused by the 
onset of XPM-effects. As Ak is increased the cascaded quadratic nonlinearity is reduced, so in 



This is a typical experimental situation, because the optimal compression point Zopt scales with A'eff [11], and since 
the nonlinear crystal length is a constant parameter one adjusts the intensity so Zopt coincides with the crystal length. 
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Fig. 3. Data from numerical simulations as in Fig.|T]for different values of Ak. (a) The FW 

is shown both for the full SEWA model and 
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duration A?opt = Af^p7"'^/1.76 at z 
when neglecting the Kerr XPM terms. The lines show the nonlocal time scales f^j, = Ti^ta.b' 
the characteristic Raman-like time Tr^shg = '^l^nl/^k, and the predicted Afopt from the 
scaling laws [ 1 1] as well as the predicted Atg°" when correcting for XPM effects on A'gff • tb 
as calculated using only up to second-order dispersion = 2, gray curve) is also shown. 
The right ordinate shows time normalized to the single-cycle pulse duration t^c = 2.0 fs. (b) 
The SHG and Kerr soliton numbers required to have A'gff = 8 fixed, achieved by adjusting 
/in. The corrected effective soliton number due to XPM effects A^^™ is also shown. 



order to keep A^eff = 8 the input intensity /in must be increased, see Fig. |3lb). Eventually the 
required intensity diverges because Ak approaches the so-called upper limit of the compression 
window Afct, beyond which A^eff < 1 always [10, 11]. As /in becomes large so does the input 
beam fluence <t>in = 2rin/in, and this makes XPM effects more pronounced; as shown in [11] 
above a critical fluence of <i>c = 33 mJ/cm^ the onset of compression in a BBO is not A^eff.c = 1 
- as one would expect from the NLS-like equation (O - but can approximately be described by 
the scaling lawA^eif.c = with AA^eff = ^'in/[*t (l +*t7<I'in)] as the delay in onset. The 

delay is caused by the XPM term creating an intensity dependent self-focusing phase-shift in 
addition to the one already created by the Kerr SPM term. The immediate consequence is that 
we can no longer expect the compression factor of fc = 33.5 predicted from Nes = 8; instead 
we must use a corrected effective soUton number N^^" = Nes — AA^eff [see Fig. Ob)], which 
for high fluences then will give a reduced compression performance. Figure [3l a) shows Afopt 
as predicted from the scaling laws for A^eff = 8 (dotted orange line), together with the corrected 
At^°" as calculated using N^^" (solid orange line). As expected for high Ak values Af™" starts 
to deviate from Afopt. Importantly, it seems to describe very accurately the compression perfor- 
mance observed numerically. An example of how the pulse looks like in this regime is shown 
in Fig.|4](Afc = 125 mm^', blue curve). The compressed pulse is longer (18 fs FWHM) and we 
also checked that it compresses later than what one would expect with A'eff = 8. These values 
correspond very well to what the reduced soliton number A^^^"^ ~ 3 predicts through the scaling 
laws. Thus, XPM strongly degrades compression when the beam fluence becomes large. This is 
also confirmed by the simulations shown with open circles in Fig. [51 a), for which XPM effects 
were turned off: only a weak degradation in compression is seen for high Ak. 

For A/: = 43 — 50 mm^ ' the limit to compression is determined by the strength of the nonlocal 



response function. Close to the transition to the nonstationary regime tf, becomes large, so the 
nonlocal response 3i+ is very broad. Initially, however, the 200 fs FWHM input pulse sees only 
a weakly nonlocal response. As the pulse compresses the nonlocal response becomes strongly 
nonlocal, whereby the NLS-like model (O reduces to a linear Schrodinger equation having a 
potential defined by the response function [19,22,24-26]. The pulse cannot be narrower than the 
width of this potential given by f^, which explains the behaviour observed for phase-mismatch 
values just above Msr. This has also been observed for spatial nonlocal solitons [19,22,24-26]. 

For < 42 mm^ ' the system is in the nonstationary regime and as Ak is reduced Afopt in- 
creases as Ak-\ The compression limit does therefore not follow th, but instead follows the 
characteristic Raman-like time shg = 2|(ii2|/Afc quite closely. Indeed some physical expla- 
nation can be extracted from this parameter, since it namely represents the pulse duration, where 
the GVM length Lgvm =A.t/\d\2 \ becomes shorter than the coherence length Lcoh = nj Nn. In- 
tuitively it seems logical that the compressed pulse duration hit a limit when the GVM length 
is equal to the coherence length: the cascaded nonlinear interaction can no longer build up the 
phase shift because the GVM will remove the FW and SH from each other before even one cas- 
caded cycle is complete. It is interesting to note that arguments similar to these were initially 
used to define the nonstationary regime [14, 18], and it was already there clear that compression 
was limited in the same way as shown here. These studies were, as mentioned before, carried 
out in the nonstationary regime as defined by the nonlocal analysis, so the results corroborate 
each other Another important result from the nonlocal analysis is that there is an oscillatory 
contribution in the nonstationary regime p (t, C/i ), see Eq. ( fTSl l. When we consider a pulse that 
gets shorter and shorter, the strength of J^[f/j^](n±) becomes larger and larger [see Fig.|2d)], 
and so we start in the present case to observe oscillations with frequency i2+ = ii'^/Tin (since 
n+ in the present example is the frequency closest to origin). That means that the soliton sees 
an oscillating potential with period so we would expect that the pulse duration can 

never go below This explains very well why the pulse duration increases in the non- 

stationary regime when ISk is decreased. It should also be mentioned that in the nonstationary 
regime Tr shg = which therefore supports this explanation. Note that no data 

points below ISk ~ 10 mm^' are shown because the cascading limit breaks down there [12]. 

To better understand the difference between the stationary and nonstationary regimes, 
Fig.Ua) shows examples of compressed pulses. For AA: = 50 mm^' a6.3 fs FWHM compressed 
pulse is observed as expected from the scaling laws, while as the transition to the nonstationary 
regime is approached (AA: = 43 mm^^) the pulse compression is limited by the nonlocal poten- 
tial of strength f/^. Once inside the nonstationary regime, the pulse not only becomes compressed 
poorly, but trailing oscillations are evident. The corresponding FW and SH wavelength-spectra 
are shown in Fig.|4|b) and (c). For M = 50 mm^' both the FW and SH spectra are very flat, 
except for a spectral FW peak and corresponding spectral SH hole. As we explain below these 
peaks are actually dispersive waves. Closer to the transition {ISk — 43 mm^') the SH spectrum 
develops a peak because here is a very narrow Lorentzian. Inside the nonstationary 

regime a distinct red-shifted peak grows up in the SH spectrum, which can be explained by the 
nonlocal theory since the spectral peak sits at the frequency In turn, close to the transition 
{ISk = 41 mm^ ' ) the FW has a corresponding spectral hole at Oi + i2+ , while further from the 
transition (AA: = 30 mm^') it becomes a spectral peak. To confirm this, we show in Fig.|4{d) 
the red-shifted holes/peaks found numerically versus AA:, with an impressive agreement with 
the nonlocal theory. This FW spectral hole/peak is the main limitation to the pulse compression 
in the nonstationary regime. 

The FW spectra in Fig. |4jb) have pronounced peaks around A = 2.9 \im. To under- 
stand this we observed that the FW dispersion operator in the frequency domain D\ (H) = 
-^QT'kY^ Chan ges sign and for Ai — 1.064 /im becomes negative beyond the dotted 
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Fig. 4. Data from the numerical results in Fig. [3] for selected values of hk: (a) \Ui'\^ at 
z = Zopt versus time, (b) the corresponding FW and (c) SH wavelength spectra. Only up to 
A = 3.5 jxm is shown in (b) since this is the edge of the transparency window of BBO [21]. 
The plots for = 50 mm^ ' and 30 mm^ ' are the cuts in Fig.[T]at z = Zopt- (d) The position 
of the red-shifted spectral peaks in the nonstationary regime. The symbols are the results 
of numerical calculations while the lines are the predictions of the nonlocal theory. 



line in Fig. |4] Based on previous results from the NLSE [27-29] at this point the compressed 
soliton should therefore be phase matched to a linear dispersive wave. Such dispersive waves 
have not been observed before with cascaded quadratic nonlinearities, but their appearance 
further underlines the analogy between soliton compression with cascaded quadratic nonlin- 
earities and with cubic nonlinear media. The wavelength of the dispersive wave Aj^ seems not 
to change as Ak is varied, which is due to the FW dispersion being independent on the crystal 
angle. Instead Ajw changes strongly with the soliton frequency, but the soliton-frequency blue 
shifts observed in the simulations were quite small and similar (around 30-40 THz). The soli- 
tons are blue shifted because M > and < and the shifts explain why the dispersive waves 
are slightly red-shifted compared to the dotted line. The dispersive waves can also be noted in 
Fig.[Tlb,d). They do not emerge before the pulse is compressed, because their strength is related 
to the spectral strength of the soliton to which they are coupled [29]. 

What happens when the effective soliton order in the stationary regime is pushed to create 
single-cycle pulses? In a previous study it was found that the GVM-induced Raman-like pertur- 
bation beyond some optimal soliton order starts to dominate and makes the compressed pulse 
asymmetric, while the peak intensity drops (Fig. 1 in Ref. [6], where Ak = 16 7r/mm). These 




Fig. 5. Results of pulse compression simulations as in Fig. [5] taking AA; = 60 mm^' and 
varying N^ff. (a) and (b) show /i, opt/An ™d Afopt when using exact dispersion (mj = oo) and 
when including up to TOD (m^ = 3), as well as the same simulations without competing 
Kerr nonlinearities. The orange curves are the predicted values from the scaling laws [11], 
corrected for XPM. (c) and (d) show FW time and spectral profiles for A'eff =17. 



results are confirmed in Fig. HJa), showing the peak intensity of the compressed pulse versus 
A'eff for = 60 mm^ ' . Beyond A'eff = 8 the pulse compression deviates from the prediction 
of the scaling law, even for the simulations including only up to third order dispersion (TOD, 
i.e., dispersion order m^/ — 3), or neglecting the competing Kerr nonlinearities. However, the 
compressed pulse duration Afopt, shown in Fig.|5tb), decreases even beyond this point of max- 
imum intensity. The explanation is that A/opt = 7in//c in (b) is determined by the compression 
factor fc alone, while the intensity in (a) is /i.opt/^in = fcQc, where Qc is the compressed pulse 
quality relating the energy of the central spike to the initial pulse energy. Around A'eff =17 the 
compressed pulse is quite close to the single-cycle regime. The time profiles for this case are 
shown In Fig.|3c). The simulations without Kerr nonlinearities (so A'shg = A'eff = 17) actually 
predict single-cycle compressed pulses, while turning on the Kerr nonlinearities the compressed 
pulses are no longer single-cycle but increase to 1.5 optical cycles. Notice also that the pulses 
with Kerr nonlinearities are more asymmetric. This asymmetry is caused by the GVM-induced 
Raman-like perturbations [1st term on the RHS of Eq. ( fT7] )l as pointed out in Ref. [6]. Since 
this effect stems from the quadratic nonlinearities it must be stressed that they are also affect- 
ing the single-cycle pulses obtained without Kerr nonlinearities. However, the difference is that 
there A'shg = 17 while with Kerr nonlinearities A'shg = 22.9 must be chosen to have A'eff =17. 
Therefore the strength of the Raman-like perturbation, which scales as N^^q, is much stronger 
when including the competing Kerr nonlinearities leading to a more asymmetric pulse. We also 
note that the simulations with exact dispersion (no polynomial expansion, = °°) have fast 



trailing oscillations, which are absent for the TOD simulations. The FW spectra in Fig.|5td) 
offer an explanation: only with exact dispersion is there a spectral peak around 3 /im; these 
are phase-matched dispersive waves. With TOD the spectrum is instead smooth because the 
phase-matching condition for the dispersive wave is pushed far into the infra-red. 

Considering these results, a brief discussion on the effect on TOD in soliton compressors 
is fruitful. In fiber soliton compressors the most detrimental effect on pulse compression is 
the Raman self-scattering term, that comes from a non-instantaneous Kerr nonlinear response. 
As mentioned before, this effect can be neglected in nonlinear crystals. However, Chan and 
Liu showed that TOD can severely distort the pulse when GVD is small [30], while for larger 
GVD the TOD is less important for the pulse shape [31], and simply leads to a slowing down 
of the soliton. In this context the TOD effect observed here, namely that the phase-matching 
wavelength of the dispersive wave is shifted, is completely different. We also remark that in our 
simulations TOD is positive, and GVD is not small {k[''^ ~ 40 fs^/mm and k^2^ ~ 100 fs^/mm). 

Summing up, the GVM-induced Raman-like distortions prevents efficient compression at 
high soliton numbers (as previously found in Ref. [6]), and the competing cubic nonlinearities 
aggravates this effect. In addition, the presence of a dispersive wave prevented soliton com- 
pression to the single-cycle level. The dispersive wave disappears when only including up to 
TOD because the dispersion is no longer accurate in the region where the dispersive wave is 
observed, and this underlines the importance of including HOD in the numerics. An interest- 
ing case to study would be when a large dispersion control is possible, such as with photonic 
crystal fibers [32]. One could namely imagine that the dispersion could be engineered, so the 
phase-matching condition for the dispersive wave is pushed into the far IR-regime, allowing for 
further compression towards the single-cycle regime. 

6. Conclusions 

To summarize we have shown that the limits to compression in cascaded quadratic soliton com- 
pressors can in most cases accurately be understood from a nonlocal model, which describes 
the cascaded quadratic nonlinearity as a nonlocal Kerr-like self -phase modulation response. 

In the stationary regime, where the nonlocal response is localized, one cannot compress 
pulses beyond the width (strength) of the nonlocal response function. When increasing the ef- 
fective soliton order to potentially compress beyond single-cycle duration, the numerical sim- 
ulations indicated that competing Kerr nonlinear effects were preventing single-cycle com- 
pressed pulses: Since the quadratic soliton number must be chosen much larger than without 
Kerr nonlinearities, this increases detrimental effects such as the GVM-induced Raman-like 
perturbation found using the nonlocal theory. Additionally it was found that higher-order dis- 
persion can also prevent the observation of single-cycle compressed pulses. In particular, dis- 
persive waves phase-matched to the compressed higher-order soliton caused trailing oscillations 
on the compressed pulse, eventually impeding further compression even at higher intensities. 

In the nonstationary regime the nonlocal response function is oscillatory, and the convolution 
between the pulse and the nonlocal response function in the weakly nonlocal limit can elegantly 
be understood in the same way as in the stationary regime, except for an oscillatory contribution 
that causes trailing oscillations in the compressed pulse and severely degrades compression. 
The SH spectrum was found to be strongly red-shifted to a wavelength accurately predicted 
by the nonlocal theory. This spectral shift in turn induces a peak in the FW, which is the main 
compression limitation in the nonstationary regime. The compression Umit was found to be 
the characteristic Raman-like response time of the cascaded process shg. roughly the pulse 
duration for which the GVM and the coherence lengths become identical. 

Another compression limit is set by the material Kerr nonlinearity, which makes compression 
possible only below a critical phase mismatch parameter, and requires large soliton orders for 



successful compression. Thus, higher-order effects (XPM, higher-order dispersion and self- 
steepening) come into play and detrimental nonlocal effects are increased. The XPM effects 
were accurately accounted for by using a corrected (reduced) soliton number 

The present analysis will serve as a useful tool for further experimental progress in soliton 
compression using cascaded quadratic nonlinearities. We will now focus our attention to com- 
pression in a BBO at Xi ~ 800 nm, because GVM is much stronger than what was presented 
here. Thus we expect the nonlocal analysis to provide more insight into this case, in particular 
concerning the nonstationary regime, which is the dominating one at 800 nm. 
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